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f ^ ■ The evolutionary conditions for the dissipative continuous magnetohydrodynamic 

^— ^ ' (MHD) shocks are studied. We modify Hada's approach in the stability analysis of the MHD 

shock waves. The matching conditions between perturbed shock structure and asymptotic 
wave modes shows that all types of the MHD shocks, including the intermediate shocks, are 
evolutionary and perturbed solutions are uniquely defined. We also adopt our formalism to 
the MHD shocks in the system with resistivity without viscosity, which is often used in nu- 
merical simulation, and show that all types of shocks that are found in the system satisfy the 
evolutionary condition and perturbed solutions are uniquely defined. These results suggest 
that the intermediate shocks may appear in reality. 

p , 

i , §1. Introduction 

o 

■ The magnetohydrodynamic (MHD) Rankine-Hugoniot relations possess six shock 

solutions; fast, slow and four types of intermediate shocks, which satisfy entropy 
condition. Arguments based upon the theory of strictly hyperbolic systems of con- 
servation laws led to the so-called evolutionary conditions (Lax 1957), and all types 
of the intermediate shocks are non-evolutionary in the ideal MHD system (Landau 
\ & Lifshitz 1960; Jeffery & Taniuti 1964; Kantrowiz & Petschek 1966; Polovin & 

Demutskii 1990; and references therein). These conditions effectively ruled out the 
intermediate shocks in nature, since they have no neighboring solution corresopnding 



> 
in 



o 



\ to small perturbations. 

Contrary to this belief, series of numerical experiments in the dissipative MHD 
. system by Wu (1987, 1988, 1990) showed that at least some of the intermediate 

shocks are admissible and can be formed through nonlinear steeping from continuous 
waves. Furthermore, Chao et al. (1993) reported the detection of an interplanetary 
intermediate shock in the Voyager 1 data. 

The dissipative steady solutions that correspond to the fast shock and the slow 
shock are coplanar (i.e., the velocity field and the magnetic field are in the same 
plane everywhere). On the other hand, four types of the dissipative steady solutions 
of the intermediate shocks can have non-coplanar structure inside their thin but 
a finite thickness front, and the non-coplanar component magnetic flux is limited 
by a maximum value, which is proportional to the dissipation coefficients. These 
properties were first pointed out by Wu (1990). 

The interactions between the intermediate shocks and the Alfven waves were 
studied by many authors through nonlinear dissipative MHD simulations (Wu 1988; 
Wu & Kennel 1992; Markovskii & Skorokhodov 2000; Falle & Komissarov 2001). 
If classical evolutionary condition is effective even in the dissipative system, the in- 
termediate shocks should instantly disappear. However, their results do not show 
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such evolution. Wu (1988) and Wu & Kennel (1992) showed that the intermediate 
shock survives a finite time after the interaction with the non-linear Alfven wave 
whose transverse magnetic field is rotated. Note that such interaction makes the 
intermediate shock non-coplanar even in the outsides of the shock structure. There- 
fore, they conjectured that there exists a new class of time-dependent intermediate 
shocks, which do not obey the MHD Rankine-Hugoniot relations, since they vio- 
late the coplanarity between the upstream and the downstream, and they are the 
neighboring states of the intermediate shocks. Wu (1988) also reported that the in- 
termediate shock remains stable as a result of the interaction with the Alfven wave. 
Therefore, the intermediate shocks seem to be evolutionary, but whether the inter- 
mediate shock evolves into other shocks and waves or not depends on the nature 
of the perturbations. Falle & Komissarov (2001) showed that the lifetime of the 
intermediate shock becomes short in a noisy circumstance. The reason is as follows. 
The magnetic flux inside the structure of the intermediate shock is rotated as a 
result of the interactions with the Alfven waves from the downstream. Then the 
intermediate shock eventually brakes up, because of the existence of maximum value 
of the non-coplanar component magnetic flux that the shock can manage. Thus, 
they cautioned that the intermediate shocks tend to survive for a long time in the 
numerical simulation, because the numerical diffusion makes the maximum value of 
the non-coplanar magnetic flux large. 

The reanalysis of the evolutionary conditions for the intermediate shocks in the 
dissipative MHD system are done by Hada (1994) and Markovskii (1998a, b). Hada 
(1994) reported that there are additional wave modes, which are originated in the 
dissipations, and the number of outgoing waves from the shock front is lager than 
the case of the ideal MHD system. He concluded that the intermediate shocks are 
evolutionary, but the set of equations are under-determined due to the excessive 
dissipative modes. As a result, he introduced the minimum dissipation principle in 
order to uniquely define solution. 

Previous authors who studied the evolutionary conditions in the dissipative 
MHD system did not take account the continuous structure of the shock front. Their 
analyses are based on the linear perturbation theory of the discontinuities or weak 
solutions. However, in the dissipative system, we have to treat the unperturbed 
shocks as a continuous transition layer, and have to solve differential equations in- 
stead of the conservation laws in the ideal system. In this paper, we examine the 
evolutionary conditions for the continuous MHD shock waves. §2 provides basic 
equations that describe perturbation of the shock structure. In §3, we formulate the 
evolutionary conditions for the continuous shock waves. We adapt our formulation 
to the various MHD shocks in §4 and §5. In §6, we summarize our result and discuss 
their implication. 

§2. Basic Equations for Linear Analysis of Dissipative MHD System 

The one-dimensional, dissipative MHD equations are described as follows: 
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where the subscript x denote the x-component and the subscript f denote the y and 
z-component. We use the unit such that the factor 4tt does not appear, v and \i are 
the shear and bulk viscosity coefficients, r\ is the electric resistivity, and k is the heat 
conduction coefficient. In addition to (|2-ip . we impose V • B = 0. This leads B x to 
a constant in the one-dimensional case. 

In the following, we consider the steady shock solution of equations (|2-ip - (|2-4p 
as an unperturbed state. In this case, without loss of generality, we can choose the 
coordinate system such that the upstream is x = —oo, the downstream is x = +oo 
(i.e., v X) q > everywhere), and the unperturbed velocity and magnetic fields are in 
the x — y plane at x = ±oo (coplanarity). We also choose the shock rest frame, and 
assume that the shock structure is around x = 0. We call the steady state type 1 if 



> Cf] type 2 if Cf > v x 



> Cf, 



type 3 if q > v x > c s ; and type 4 if c s 



> v x , where 



Cf, c s and q are the fast, slow and Alfven (intermediate) speed, respectively. Steady 
shock solutions of equations (|2-ip - (|2-4p are well studied (see, e.g. Wu 1990). There 
are six types of the shock solutions, 1 — ► 2 fast shock, 3 — ► 4 slow shock and 1 — > 3, 
1 — > 4, 2 — ► 3 and 2 — > 4 intermediate shocks, where the numbers before and after 
the arrow represent the state of ahead and behind the shock front. 

Let us consider the small perturbation of the steady shock structure. We assume 
that the perturbation of the physical variable g(x, t) takes the following form: 



g(x,t) = g (x) + 8g(x)e 



-i u>t 



(2-5) 



where the subscript denotes the unperturbed variable. Linearizing the equation 
(|2Tp . we obtain the perturbed shock equations 
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—Sv z = — (v xfi v zfi Sp + po v zfi Sv x + po v xfi 8v z - B xfi SB Z ) 



-i to (p Sv x + v X:0 Sp) , (2-7) 
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+-B )Sv x - (5v x B xfi + 5v y B yfi + Sv z B z>0 + v y SB yfi 
+v z 6B z>0 + )B x } -iujSe. (2-12) 



Equation (|2-6p is the first-order ordinary differential equation with respect to 
Sp, and equations (|2-7j) - (|2-12j) are the second-order differential equation with respect 
to 5v x , Sv y , 5v z , SB y , SB Z , and Sp respectively. Therefore, we have totally the 13th 
order differential equations. We denote the order of the perturbed equations as 
TV" = 13. If we choose the coplanar intermediate shock solution (v z o = -B 2 ,o = 
everywhere) as an unperturbed state, equations (|2-6p - (j2-12p are separated into the 
two sets of equations. The first set is equations (|2-6p - (|2-8|) . (|2-10|) . and (|2-12|) . which 
relate the x and y- components of the perturbations, and they are totally the 9th 
order differential equations (N xy = 9). The other set is equations (|2-9p and (|2Tip . 
which relate the ^-components of the perturbations, and they are totally the 4th 
order differential equations (N z = 4). 

§3. Evolutionary Conditions for the Continuous Shock Waves 

In this section, we formulate the evolutionary conditions for the continuous MHD 
shock waves. We divide the space into three regions. Two of them are the regions 
far away from the shock front, in which the unperturbed physical variables can be 
regarded as constants. We call such regions in the upstream and the downstream 
as region IA and T>, respectively. The region between the regions U and T> is the 
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transition region where the unperturbed physical variables change continuously. We 
call it region T. The situation is schematically illustrated in Figure [TJ 
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Fig. 1. Schematic of region partition. 



As shown in the previous section, perturbed shock equations are totally iVth 
order ordinary differential equations. Thus, in order to determine a solution of these 
equations in the region T, we need N boundary conditions at the edges of the region 
T. In other words, in the region T, there are N degree of freedom in order to 
determine a solution of these equations. 

In the regions IA and V, unperturbed physical variables are regarded as con- 
stants, and we can obtain the asymptotic solutions of the perturbed shock equations. 
Omitting the spatial derivative of the zeroth order variable, and Fourier transform- 
ing equations (|2-6[) - (j2- 12|) in space id/dx ~ ik), one can obtain the characteristic 
equation for the asymptotic waves. The determinant of the characteristic matrix 
provides the characteristic equation. 
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V' = k 2 {2,p-Au)/p-2,iQ 
R = r) -iQ/k 2 , 



(3-9) 
(3-10) 
(3-11) 
(3-12) 



K = (7- l) K /p-iQ/k 2 , 
K' = k 2 {-/-l)K/p-ij[2, 



where the zoroth order variables are evaluated in the region Li and T>, and we use the 
fact that the non-coplanar component of the unperturbed variables vanish {B z $ = 
v z,o = 0) in the asymptotic regions owing to the coplanarity. The solutions of the 
characteristic equation Df se = for uj as a function of k provide the dispersion 
relations of the fast, slow, and entropy waves in the uniform dissipative medium, 
and the solutions of Di = provide the dispersion relation of the Alfven waves. 

Definition of Mode As shown below, we use the asymptotic solutions in 
regions U and P as a boundary conditions of the differential equations (|2-6j) - ([2-12p . 
These conditions are expressed as the superpositions of the independent asymptotic 
solutions. The degree of freedom of the spatial behavior of the perturbation for given 
uj gives the number of independent asymptotic solutions. Thus, they are obtained 
by solving characteristic equation (|3-1|) for k as a function of given uj. In this paper, 
we call it "mode" or "asymptotic mode". Note that it does not mean solution of 
characteristic equation for uj as a function of k. 

Hada (1994) studied the solutions of D{ = in the limit of small dissipation 
coefficients. He found that in addition to the solutions, which correspond to the 
Alfven modes, there are additional dissipative modes. Details of the solutions are 
discussed in the next section. 

Let us consider evolutionary condition. If we steadily throw a small amplitude 
incident (ingoing) wave whose frequency is u> toward the shock from the region Li or 
P, then the shock front is perturbed and emit the waves with the same frequency. 
Let m (= mu + mx>) denote the number of resulting asymptotic modes in the region 
U {mu) an d T> {mx>) without the incident wave, i.e. the number of modes that can 
be emitted or raised at the shock front. In the regions hi and V, the solutions of 
perturbed shock equations should be expressed as the superpositions of the m asymp- 
totic modes and one incident wave. These asymptotic solutions are determined by 
m parameters, i.e. amplitudes of the m asymptotic modes. Note that the amplitude 
of the incident wave is determined. Thus, in the case that 



we can obtain the solution of the perturbed shock equations ()2-6j) - f)2- 12|) in the region 
T, which is smoothly connected to the asymptotic solutions at the edges of the region 
T, for arbitrary incident wave by choosing m amplitudes of the asymptotic modes. 

If m is less than N, we cannot obtain solution. If m is greater than N, we can 
obtain solution, but we need additional conditions or constraints in order to uniquely 
define the solution. Therefore, the condition (|3T3p corresponds to the evolutionary 
condition in the dissipative system. 



N = m 



(3-13) 



Evolutionary Conditions for MHD Shock Waves 



7 



If unperturbed shock structure is coplanar, the evolutionary condition (13T3|) can 
be divided into the following two conditions 

N xy = m xy , (3-14) 

N z =m z , (3-15) 

where m xy and m z are the numbers of the asymptotic modes of Df se = and D{ = 0, 
respectively. We call the condition (|3- 14j) evolutionary condition for the x and y- 
components, and condition (|3-15j) evolutionary condition for the ^-components. Of 
course, evolutionary shock must satisfy both conditions. 

§4. Evolutionary Condition for the z-coraponents 

In this section we show that, contrary to the ideal system, all types of inter- 
mediate shocks satisfy evolutionary condition for the z-components N z = m z . The 
solutions of Di = for k under the assumption of small dissipation coefficients, 
derived by Hada (1994), are 

k^ = —^— + 0( V ,u), (4-1) 

V x ,0 ± Ci 

k (±d) = J^((V + *>) v*,o ±{(V + vfvlo -4171/ «o - <$ )}*) + O( V , u \4-2) 

The (+) and (— ) modes are the usual Alfven waves, which propagate parallel and 
anti-parallel to the a;- axis in the fluid rest frame. The (±<i) modes, with no coun- 
terparts in the ideal system, are so called the "dissipative modes". The dissipative 
modes do not propagate, because k^ d ^ do not have real part in their primary terms, 
and they are localized around the shock front. We exclude diverging dissipative 
modes as unphysical asymptotic solutions. Thus, they are physical if 

Im[k {±d) ] > in the region U, or (4-3) 

Im[fc (±d) ] < in the region V. (4-4) 

If the flow speed is super- Alfvenic in the shock rest frame, then from equation (|4-'2|) . 
the imaginary part of k^ d > are less than zero. Therefore there are two dissipative 
modes in the region U, whereas no dissipative mode in the region D. On the other 
hand, if the flow speed is sub- Alfvenic, the imaginary part of k^ +d ^ is greater than 
zero and that of k^~ d ^ is less than zero, thus there is one dissipative mode in the 
region U and T>. 

Let us consider the evolutionary condition for the z-components. In the case of 
the fast shock, the flow speeds both in the upstream and the downstream are super- 
Alfvenic. Therefore, the fcW modes propagate to the downstream in the region T>, 
and the k^ d ^ modes are in the region IA. In the same way, the flow speeds of the 
slow shock are sub- Alfvenic in the both sides. There are and k^ +d ^ modes in the 
region IA, and k^ and k^~ d ^ modes in the region T>. In the case of the intermediate 
shocks, the upstream flow speed is super-Alfvenic and the downstream flow speed 
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is sub-Alfvenic. There are k^ d ^ modes in the region U, and k^ and k^~ d ' modes 
in the region T>. The situations are schematically illustrated in Figure [2j In all 
cases, the MHD shocks satisfy the evolutionary condition for the z-components, i.e. 
N z = m z = 4. 

downstream 

V x > Ci 

(-) Alfven 
(+) Alfven 

V x < Ci 

(+d) dissipative 
(+) Alfven 

V x < Ci 

(+d) dissipative 
(+) Alfven 

Fig. 2. Illustration of the outgoing ordinary MHD modes and the localized dissipative modes, where 
the bars represent the region T, the arrows represent the propagation of the asymptotic modes, 
and the triangles represent the localized dissipative modes. The (±) sign before the Alfven 
modes denote the propagation direction in the fluid rest frame, and c; = (-B^o/Po) 1 ^ 2 is the 
Alfven (intermediate) speed. 



§5. Evolutionary Conditions in the Resistive but Inviscid MHD System 

Let us consider the case where there is only resistivity in the dissipation. The 
one-dimensional basic equation of the weakly ionized gas, a main ingredient of the 
interstellar medium, can be written as a resistive but inviscid MHD equation un- 
der the strong coupling (one- fluid) approximation. Furthermore, Wu (1987, 1990) 
showed by using one-dimensional simulation in this system that the intermediate 
shocks are formed through nonlinear steeping from simple waves. Thus, it is mean- 
ingful to analyze the evolutionary condition in this system. We apply our formalism 
developed in §3 to the resistive but inviscid MHD shocks. 

In this system the perturbed basic equations ()2-6j) - (|2- 12[) become first order 
ordinary differential equations with respect to dp, 5p, 5v, and second order one with 
respect to 5Bt- Thus, we have totally 9-th order (N = 9) differential equations 
(N xy = 6, N z = 3 in the case of the coplanar shocks). The number of asymptotic 
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modes is determined by the solutions of the resistive version of equation (13-111 . 

D fse = Q {c 2 c 2 a k 4 - in c 2 a Q k 4 - (c\ + c 2 a )Q 2 k 2 + z 77 Q 3 k 2 + Q 4 } (5-1) 
Di = Q(Q + ink 2 ) - c 2 ifl k 2 . (5-2) 

Equation ()5 2|) . which contains three asymptotic modes, was studied by Hada (1994). 
Two of them corresponds to two Alfven waves, and third one corresponds to the 
dissipative mode. The dissipative mode solution under the assumption of small 
resistivity is 

k (di) =l c "~<o +Q( o o) (5.3) 
V v x ,o 

which is physical in the upstream for the super- Alfvenic state (state 1 and 2), and 
physical in the downstream for the sub- Alfvenic state (state 3 and 4). 

Equation (15-ip contains six asymptotic modes. Five of which correspond to two 
fast waves, two slow waves and one entropy wave, and the last one corresponds to the 
dissipative mode. The dissipative solution under the assumption of small resistivity 
is 

k(d2) = . c 2 a c 2 -(c 2 A + c 2 a )vl + vl + ^ 

In the case of the MHD system with resistivity without viscosity, steady shock struc- 
ture whose flow velocity is supersonic (v x > c a ) ahead of the shock front and subsonic 
(v x < c a ) behind the shock front is impossible (see, Wu 1990). Thus, it is possible 
that the shock structure for 1 — ► 2, 1 — ► 3, and 2 — > 3 whose flow speeds are super- 
sonic everywhere in the shock rest frame, and 2 — > 3, 2 — > 4, and 3 — > 4 whose flow 
speeds are subsonic everywhere in the shock rest frame. The denominator of the 
righthand side of equation (|5-4|) is less than zero in the supersonic state, and larger 
than zero in the subsonic state. The numerator is larger than zero in the state 1 
(super-fast) and 4 (sub-slow), and less than zero in the state 2 and 3 (sub-fast and 
super-slow). Therefore, in the supersonic case, the dissipative (d2) mode is physical 
in the region U for the state 1, and in the region T> for the state 2 and 3. In the 
subsonic case, it is physical in the IA for the state 2 and 3, and in the T> for the state 
4. 

The number of asymptotic modes in asymptotic regions is m = 9 {m xy = 
6, m z = 3 in the case of non-coplanar shock) for all shocks. We illustrate the situa- 
tions for the shocks whose upstream velocity is supersonic in Figure [31 and for the 
shocks whose upstream velocity is subsonic in Figure HI Therefore, all the shocks, 
which are allowed in the MHD system with resistivity without viscosity, satisfy the 
evolutionary condition N = m = 9 (or N xy = m xy = 6, N z = m z = 3). 

§6. Summary and Discussion 



We have studied evolutionary conditions for the continuous MHD shock waves 
in the dissipative system. We have shown that all types of MHD shocks, even the 
intermediate shocks, satisfy the evolutionary condition for the z-components, which 
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upstream -pg downstream 

(-) Alfven 
(+) Alfven 
(d2) dissipative 
(-) slow 
entropy 
(+) slow 
(+) fast 



(dl) dissipative 
(+) Alfven 
(d2) dissipative 
(-) slow 
entropy 
(+) slow 
(+) fast 



(dl) dissipative 
(+) Alfven 
(d2) dissipative 

(-) slow 
entropy 
(+) slow 
(+) fast 

Fig. 3. Illustration of the outgoing ordinary MHD modes and the localized dissipative modes 
for the shocks whose upstream velocity is supersonic, where the bars represent the region T, 
the arrows represent the propagation of the asymptotic modes, and the triangles represent 
the localized dissipative modes. The (±) sign before the ordinary MHD modes denote the 
propagation direction in the fluid rest frame. 

is not satisfied in the ideal MHD system. Especially in the resistive system, all types 
of MHD shocks, which are allowed in the resistive system, satisfy the evolutionary 
condition. Therefore, the intermediate shocks have neighboring solution, and they 
can survive interactions with arbitrary small amplitude waves. 

The difference of the evolutionary conditions between the ideal and the dissi- 
pative system is that, in the dissipative system, there are dissipative perturbation 
modes, which do not propagate and are localized around the shock front, and they 
provide degree of freedom, which are lost in the ideal system. This difference is 
first pointed out by Hada (1994). Hada (1994) tried to connect the perturbation by 
using the Rankine-Hugoniot relations. As a result, set of equations become under- 
determined due to the excessive dissipative modes. In order to uniquely define so- 
lution, Hada (1994) introduced the minimum dissipation principle. In this paper, 
however, we have shown that perturbations in the asymptotic regions (far away from 
the shock front) can be uniquely connected by solving differential equations, instead 
of the conservation laws in the ideal or hyperbolic system, because the shock is not 
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(dl) dissipative 
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(-) fast 
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downstream 

(dl) dissipative 
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(-) slow 
entropy 
(+) slow 
(+) fast 



(dl) dissipative 
(+) Alfven 
(d2) dissipative 
entropy 
(+) slow 
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(dl) dissipative 
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(d2) dissipative 
entropy 
(+) slow 
(+) fast 



Fig. 4. Illustration of the outgoing ordinary MHD modes and the localized dissipative modes for the 
shocks whose upstream velocity is subsonic, where the bars represent the region T, the arrows 
represent the propagation of the asymptotic modes, and the triangles represent the localized 
dissipative modes. The (±) sign before the ordinary MHD modes denote the propagation 
direction in the fluid rest frame. 



a weak solution but a continuous solution in the dissipative system. Note that the 
solution of the differential equations also satisfies the usual conservation laws in the 
asymptotic regions. 

However, the idea of the minimum dissipation principle and use of the perturbed 
Rankine-Hugoniot relations may be important, when we analyze the linear stability 
of the intermediate shocks. In order to analyze the linear stability by using our 
formulation, we have to solve the differential equation. This seems difficult, because 
the equations are stiff, and there need huge amount of calculations to search param- 
eter space. If the minimum dissipation principle can sort out important dissipative 
modes from other dissipative modes in order to the solution be unique, the analysis 
is easily done, because the solution is obtained analytically. Further study is needed 
to discover such technique. 

Our approach is similar to that of Wu (1988) and Wu &; Kennel (1992). They 
studied the interaction between the intermediate shock and the non-linear Alfven 
wave whose transverse magnetic field is rotated. They showed by using non- linear 
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simulations that there exists a new class of time-dependent intermediate shocks, since 
they violate the coplanarity even without the shock structure (i.e., they do not obey 
the MHD Rankine-Hugoniot relations between the upstream and the downstream), 
which possibly be a neighboring state of the intermediate shocks. On the other hand, 
our approach treat the interaction between the intermediate shock and small ampli- 
tude (linear) MHD wave, and showed that the intermediate shocks have neighboring 
solutions which describe the time evolutions of the perturbed intermediate shocks in 
the linear regime. 

The intermediate shock which has maximum non-coplanar magnetic flux inside 
the shock structure break up into other shocks and waves as a result of the interaction 
with the Alfven wave (Markovskii & Skorokhodov 2000 and Falle & Komissarov 
2001). Such an unstable phenomenon would be analyzed possibly in terms of the 
linear stability analysis. Our formulation is also expected to be used for such analysis. 

The existence of the intermediate shocks has been argued mainly in the inter- 
planetary system (Chao et al. 1993; Chao 1995). However, in the partially ionized 
interstellar medium, it is known that there are C-type MHD shocks, which have 
broad shock structure owing to friction between neutral and ionized gases. Because 
of the lifetime of the intermediate shock seems to be proportional to the width of 
the front, the intermediate shocks can be long-lived in the interstellar medium. The 
quantitative discussion on this issue remains to be done. 
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